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The possible existence of maximal or near-maximal lepton mixing con- 
stitutes an intriguing challenge for fundamental theories of flavour. We study 
the phenomenological consequences of maximal and near-maximal mixing of 
the electron neutrino with other (x=tau and/or muon) neutrinos. We describe 
the deviations from maximal mixing in terms of a parameter e = 1 — 2 sin^ O^-x 
and quantify the present experimental status for |e| < 0.3. We show that 
both probabilities and observables depend on e quadratically when effects are 
due to vacuum oscillations and they depend on e linearly if matter effects 
dominate. The most important information on i/g -mixing comes from so- 
lar neutrino experiments. We find that the global analysis of solar neutrino 
data allows maximal mixing with confidence level better than 99% for 10~^ 
eV^ < Am^ < 2 X IQ-'^ eV^. In the mass ranges Am^ > 1.5 X 10"^ eV^ 
and 4 x 10"^° eV^ < Am^ < 2 x 10~^ eV^ the fuh interval |e| < 0.3 is al- 
lowed within ~ 4cr (99.995 % CL). We suggest ways to measure e in future 
experiments. The observable that is most sensitive to e is the rate [NC]/[CC] 
in combination with the Day-Night asymmetry in the SNO detector. With 
theoretical and statistical uncertainties, the expected accuracy after 5 years 
is Ae ~ 0.07. We also discuss the effects of maximal and near-maximal 
fg-mixing in atmospheric neutrinos, supernova neutrinos, and neutrinoless 
double beta decay. 
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I. INTRODUCTION 



The data from both atmospheric and solar neutrino experiments have given a rather con- 
vincing evidence for non-zero neutrino masses and mixing. Furthermore, some features of 
the neutrino flavour parameters are surprising and quite different from the quark flavour pa- 
rameters. In particular, one or two of the three mixing angles in the MNS (Maki-Nakagawa- 
Sakata) mixing matrix for leptons [1] are large. Specifically, the simplest interpretation of 
the atmospheric neutrino measurements gives [2,3] 

sin^ 29^r ~ 0.8 - 1, Am^2 ~ (2 - 5) x 10~^ eV^. (1) 

There exist several solutions of the solar neutrino problem involving oscillations of electron 
neutrinos into some combination, u^, of /x— and r— neutrinos with large mixing angle [4-10] 
with parameters in the range 

sin^ 29ex - 0.7 - 1, Amj^ - (0.2 - 4) x 10"^ eV^ or (0.05 - 20) x 10"^ eV^ (2) 

The mixing angles involved in the explanation of the solar and atmospheric neutrino 
data are not just order one. They are actually near-maximal, that is, sin^2^ close to 1. If 
indeed one of the mixing angles is near-maximal, it would provide a strong support to the 
idea that the corresponding neutrinos are in a pseudo-Dirac state. Such a scenario would 
have very interesting implications for theoretical flavour models. These implications have 
been recently studied in Ref. [11]. A precise knowledge of the mixing and, in particular 
lower and/or upper bounds on small deviations from maximal mixing provides an excellent 
probe of the related flavour physics. It is the purpose of this work to study the experimental 
ways in which the region of near-maximal mixing can be probed. 

It is important to notice that there exists also a viable solution of the solar neutrino 
problem that does not involve large mixing, the so called small mixing angle solution (SMA). 
Clearly, identiflcation of the SMA solution will immediately exclude the possibility discussed 
in this paper. Also discovery of the sterile neutrinos will require a change of the whole 
picture. Here we consider only three light active neutrinos and take the 3x3 MNS matrix 
to be unitary. 

Our main interest lies in the study of near-maximal mixing involving Ug. We deflne a 
small parameter e which describes the deviation from maximal mixing as: 

sin2^ex = ^(l-e), (3) 

Our convention is such that Am|i = m| — > 0. Then e > (< 0) corresponds to the 
case that the lighter neutrino mass eigenstate, ui, has a larger (smaller) component of Ve 
and the heavier one, U2, has a larger (smaller) component of u^. 

Which value of deviation from maximal mixing is expected? In the case of theoretical 
models where the pseudo-Dirac structure is naturally induced, one expects that the deviation 
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from maximal mixing is suppressed by a small parameter that is related to an approximate 
horizontal symmetry. If the same symmetry is also responsible for the smallness and hi- 
erarchy of the quark sector parameters, then it is quite plausible that |e| < C(A), where 
A = 0.22 is the Cabibbo angle in the Wolfenstein parametrization. In various flavor models, 
the deviation from maximal mixing is related to other physical parameters. For example, in 
a large class of models, we have [12] 

|e| > 2me/m^ ^ 0.01. (4) 

The U{1) X f/(l) models described in refs. [13,11] give |e| ~ {Aml2/ ^'mlsV^^ , while models 
of Lg — — Lt symmetry [14,11] give |e| ~ Am^2/^^23- Quite generally we have |e| > 
Am^2/'^^23' which could be more restrictive than Eq. (4) if the solution to the solar neutrino 
problem lies at the upper end of the Am^2 range given in Eq. (2). In a large class of models 
we also have |e| of the same order of magnitude as jf/esl, the mixing of Ue in the third mass 
eigenstate, which can again be tested in the future [11]. 

Large values of e, |e| > 0.3, will testify against at least the simplest versions of these 
theories. Therefore we consider both positive and negative values of e in the range 

|e|^0.3 (5) 

which corresponds to sin^ 26ex ^ 0.91. As concerns the mass-squared difference we cover the 
whole range below the reactor bound: 

10-^^ eV^ < Am^i < lO'^ eVl (6) 

Most of our discussion takes place in an effective two neutrino generation framework. 
This is justified if U^^ is zero or sufficiently small. We quantify this condition and consider 
the effect of a non-zero U^.^ in the Sect. VI. We find there that, for matter oscillations, the 
leading corrections to the case of maximal mixing are of C(e, If/esP), so that a reduction to 
a two neutrino analysis is justified for \Uez\'^ <C e. For vacuum oscillations, the corrections 
are of 0(e^, |f/e3p), and a two neutrino analysis is valid for \Ue^\ <^ e. 

On the experimental side, wc note that if Am^^^^ > 2 x 10^^ eV^, then the CHOOZ ex- 
periment limit [16] requires |^7e3p < 0.05 [17,6]. For such values of |?7e3p, a two generation 
analysis of matter (vacuum) effects is always valid for |e| > 0.1(0.3). On the theoretical 
side, in a large class of flavour models, \Uez\ ~ ]e| [11]. In such a framework, a two genera- 
tion analysis of matter effects is a good approximation while vacuum oscillations should be 
considered in the three generation framework. 

In the limit \Ue.z\ = (which reduces the problem to a two neutrino framework), the 
deviation from maximal mixing can be determined as in (3). In addition the mixing of 
and Uj- can also be maximal, as favored by the atmospheric neutrino data. In this case we 
have the mixing structure 

\Ues\' = 0, \Ue,\' = \U,,\^ = |t;^3|2^|t;^3|2^ 1^ (7) 

which corresponds to the so called bi- maximal mixing scheme [15]. The analysis for the 
solar neutrino phenomenology is independent of [/^3 and therefore the results discussed 
in Sees. Ill, IV, V, VII B and VII C apply also to the bi-maximal mixing case. Only in 
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Sec. VII A where we discuss atmospheric neutrinos, does Ufj^s play a role, and there we take 
it to be large and possibly maximal. 

The plan of this paper is as follows. In Sec. II we give some basic physics considerations 
and useful expressions for the survival probability and for various observables. In Sec. Ill we 
describe the present status of maximal mixing from solar neutrino experiments. The results 
of a global fit to all available solar neutrino data are given in Sec. Ill A. The dependence 
of these results on various aspects of the analysis is described in Sec. IIIB. In Sec. IV we 
study the predictions that follow from near-maximal mixing for individual, existing mea- 
surements: total rates. Argon production rate. Germanium production rate, the Day-Night 
asymmetry in elastic scattering events, the zenith angle distribution of elastic scattering 
events, and the shape of the recoil energy spectrum. In Sec. V we suggest tests of maxi- 
mal mixing from future experiments. We describe the conditions for having unambiguous 
tests in Sec. VA. Then we examine individual experiments: GNO and Super-Kamiokande, 
SNO, Borexino and low energy solar neutrino experiments. The effect of extending to three 
neutrino scenario is commented in Sec. VI. In Sec. VII we discuss the effect of maximal 
mixing in atmospheric neutrinos, supernova neutrinos and neutrinoless double beta decay. 
We present our conclusions in Sec. VIII. 

II. PHYSICS AT NEAR-MAXIMAL MIXING 
A. The survival probability for solar neutrinos 

In this subsection we present general expressions for the survival probability of solar 
electron neutrino in a two generation framework valid in the full range of Am^ which we 
use in our numerical calculations. 

The survival amplitude for a solar Ue neutrino of energy E a.t a. detector in the Earth 
can be written as: 

Aee = E 4. eM-im^L - r)/2E] . (8) 
1=1 

Here A^^ is the amplitude of the transition z/g — > z/j (z/j is the i-mass eigenstate) from the 
production point to the Sun surface, Af^ is the amplitude of the transition z/j z/g from the 
Earth surface to the detector, and the propagation in vacuum from the Sun to the surface 
of the Earth is given by the exponential. L is the distance between the center of the Sun 
and the surface of the Earth, and r is the distance between the neutrino production point 
and the surface of the Sun. The corresponding survival probability Pee is then given by: 

Pee PlPle + P2P2e + 2^fP^^P^^ COS ^ 

= Pi + (1 - 2Pi)P2e + 2^Pl(l - Pl)P2e(l - P2e) COS^. (9) 

Here Pi = |v4fjp is the probability that the solar neutrinos reach the surface of the Sun as 
and we use Pi -I-P2 = 1, while Pje = \AfJ'^ is the probability of Ui arriving at the surface 
of the Earth to be detected as a i^e, and we use Pie + P2e = 1- The phase ^ is given by 
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(10) 



where 5 contains the phases due to propagation in the Sun and in the Earth and can be 
safely neglected since it is always much smaller than the preceeding term [18,19]. 
From Eq. (9) one can recover more familiar expressions for Pee' 

(1) For Am? / E < 5 x 10"" eV, the matter effect supresses flavour transitions in both the 
Sun and the Earth. Consequently, the probabilities Pi and are simply the projections 
of the z/g state onto the mass eigenstates: Pi ~ cos^ ^, P-ie — sin^ Q. In this case we are left 
with the standard vacuum oscillation formula: 

PJ^^l- sm^2esm'^(Am^(L - r)/AE) (11) 

which describes the oscillations on the way from the surface of the Sun to the surface of the 
Earth. The probability is symmetric under ^ <-> | — ^. 

(2) For Arm? /E > 10~^^ eV, the last term in Eq. (9) vanishes and we recover the incoher- 
ent MSW survival probability. For Air? /E ~ 10~^^ — 10~^^ eV^, this term is zero because 
Ve adiabatically converts to V2 and Pi = 0. For Air? /E ^ 10~^^ eV^, both Pi and P2 are 
nonzero and the term vanishes due to averaging of cos^. 

(3) In the intermediate range, 5 x 10-^^ < An? IE < 10-^^ eV, adiabaticity is violated 
and the cos^ coherent term should be taken into account. The result is similar to vac- 
uum oscillations but with small matter corrections. We define this case as quasi-vacuum 
osciUations [18-20]. 

The results presented in the following sections have been obtained using the general 
expression for the survival probability in Eq. (9) with Pi and found by numerically solving 
the evolution equation in the Sun and the Earth matter. For P, we use the electron number 
density of BP2000 model [21]. For Pj^ we integrate numerically the evolution equation in 
the Earth matter using the Earth density profile given in the Preliminary Reference Earth 
Model (PREM) [22]. 



B. The mixing angle in matter 

While, as explained above, our results are obtained by a numerical calculation, it is 
useful to find approximate analytical expressions for the neutrino survival probability and 
for various observables. This is done in this and in the following subsections. The analytical 
expressions help us to qualitatively understand the behaviour of the different observables, 
particularly for the case of near-maximal mixing. 

The probabilities and observables depend on the mixing angle in matter 9m via cos 29^, 
which enters the probability of the adiabatic conversion and via sin^ 29jn which determines, 
e.g., the depth of oscillations in a uniform medium. 

We can write cos 26 ^ in terms of the neutrino oscillation parameters and the electron 
density in medium as: 

— 1+1] cos 29 , . 

(1 — 2r] cos 29 + rj^yi^ 

Here 
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is the ratio between the refraction length, Iq, and the neutrino oscillation length in vacuum, 

27rmjv , 4:TrE 

Here p is the matter density and is the number of electrons per nucleon. 
Around maximal mixing cos 29^ can be expanded as 

-^(l-^.). (15) 

In the limit of weak matter effects, 77 ^ 1, and in the matter dominance case, 77 <^ 1, we 
get: 

cos2^^ = |^_^ ;;<<}■ (16) 

The dependence of cos 29^^ on e is smooth. It is stronger for 77 ^ 1 and highly suppressed 
for <S 1. For precisely maximal mixing we have cos2^j„ = —l/^/l^^rf, which decreases 
from zero in vacuum to —1 in the matter dominance case. 

The expression for sin^ 29m for near maximal mixing is given by 

-^29.^^(l + ^e), (17) 

which leads to 

sin^2^„=|M^, , . (18) 

\ 7/2(1 + 277e) 77 < 1 ^ ^ 

In both cases the e corrections to the value at maximal mixing are strongly suppressed. In 
vacuum, the correction is quadratic in e: sin^ 26* = 1 — e^. 



C. Survival probability 

We first consider the survival probability for electron neutrinos without the regeneration 
effect in the Earth. It describes the flux arriving at the Earth during the day. In daytime, 
P2e = sin^ 9. Consequently, Eq. (9) gives, in the region where the oscillating term in Eq. (9) 
is absent. 

The neutrino evolution in the Sun described by the probability Pi can be approximated by 
the well known expression 
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A = ^+Q-Pc)cos2^5. (20) 

Here Og is the matter mixing angle at the production point: 

COS26I5 = cos26'^(775), r]s = r]{psYes), (21) 

where ps and Y^s are, respectively, the matter density and the number of electrons per 
nucleon at the production point. Eq. (20) is assumed to be averaged over the production 
region in the Sun. Pc is the jump probability which describes the violation of adiabaticity. 
For an exponential density profile it takes the following form [23,24]: 

-Tsin^e _ -7 g-(7/2)(l-e) _ -7 

Pc = = , (22) 

l_e-7 i_e-T ' ^ ^ 

where 7 is the ratio of the density scale height Ip and the neutrino oscillation length: 

The length scale ro = i?0/lO.54 is related to the exponential approximation to the solar 
density profile, p = po exp(— r/ro). 

Notice that, originally, Eq. (22) was derived for a mixing angle 6 < j where resonant 
enhancement is possible. However both Eq. (12) and Eq. (22) can be analytically continued 
into the second octant, 9 > 7r/4, and used to compute the corresponding survival probability 
for e < [10,5-7]. 

Inserting Pi of Eq. (20) into Eq. (19), we get 

PD = ^ + e(l-Pc)cos29s. (24) 



2 V2 

Let us study the properties of Pp. In Fig. 1 we plot Pd as a function of Am^/AE for 
different values of the deviation from maximal mixing. We show in the figure the results 
obtained by the numerical calculation as well as from the corresponding analytical approxi- 
mation (22) for exponential density profile. We learn the following points from Eq. (24) and 
Fig. 1: 

1) For maximal mixing (e = 0), Pd — 1/2 independently of the adiabaticity violation 
(encoded in Pc), matter effects, etc.. 

2) For near-maximal mixing (e 7^ 0), solar matter effects lead to an energy dependent 
probability in the MSW region 4 x 10"^^ < Am'^/E < 10"^° eV. 

3) For Arri^/E > 10"^^ eV, the adiabaticity condition, 7 > 1, is fulfilled (see Eq. (23)). 
Consequently Pc — ^ and Eq. (24) gives 

= - (1 + e cos 2^5) {Am^/E > 10"^^ eV). (25) 

From Fig. 1 we see that the analytical expression gives a good description of the propagation 
in the Sun for this region. 

4) For 3> 1 (weak matter effect), we have cos 29s — ^ (see Eq. (16)) and the probability 
(25) reduces to the vacuum oscillation probability: 
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(26) 



5) For ?75 -C 1 (the resonance layer is far enough from the production point which is 
close to the center of the sun), we have cos 2^5 = — 1 and 



This is the region of strong adiabatic conversion, when v^. produced in the center of the Sun 
practically coincides with the matter eigenstate at the production point and during its 
propagation inside the Sun. Therefore, it emerges from the Sun and reaches Earth as z/2- 

6) For ^w? / E < 10~^^ eV, adiabaticity is violated. We see from Fig. 1 that the analytical 
results obtained for an exponential density profile differ from the results of the numerical 
calculations. In particular, the analytical result shows a "slower" transition to the vacuum 
oscillation regime or, in other words, it overestimates the size of the matter effects in the 
quasi-vacuum oscillation region. The same value of the survival probability appears for 
about two times smaller lS.w? jE. Similar conclusions have been drawn in refs. [19,20] . 

7) For e 7^ 0, the effects of the adiabatic edge situated at AmV-E = (lO^^^ - 10"^°) 
eV, and of the non-adiabatic edge situated at Am^/£^ = (10~^^ — 3 x 10"''^^) eV, become 
important. 

8) For !\vp?jE <^ 10^^^ eV, as we noted in Sec. II A, the effect is reduced to the 
vacuum oscillation between the surfaces of the Sun and the Earth and the average survival 
probability, Pd = 1 — | sin^ 29, (shown in Fig. 1) coincides with that in Eq. (26). However, 
in this region averaging of oscillations does not occur and for the survival probability we 
should use 



where = Am^(L — r)/4:E is the oscillation phase which does not depend on e. In principle 
matter effects strongly suppress the oscillations inside the Sun and the Earth. However, the 
modification of the observables is small, since the size of the Sun (and the Earth) is much 
smaller than the oscillation length in vacuum. The quadratic dependence of the probability 
on e is again restored. Moreover, time variations of signals and the distortion of the energy 
spectrum originate from the dependence of the phase on L and E. Therefore according 
to (28), both the variations and the distortion are proportional to (1 — e^). That is, the 
dependence of all observables on near maximal mixing is very weak. 

According to Eq. (27) and Fig. 1, in the MSW region (more precisely at the bottom of 
the suppression pit), the survival probability depends on e linearly. For Am^/E below and 
above the MSW region the probability converges to the vacuum oscillation probability. The 
deviation of the latter from the probability at maximal mixing depends on e quadratically 
and the dependence on the sign of e disappears. From this we infer that the sensitivity of 
Pd and consequently of observables to e is much higher in the MSW region. 

Notice that there are two transition regions (between MSW and the vacuum oscillation 
regions) where the effects are mainly due to vacuum oscillations with small matter cor- 
rections. We call these regions 'Quasi- Vacuum Oscillation regions': we denote by QVOl 
(QVOs) the one with Am^ larger (smaller) than in the MSW region. In the QVOl and 
QVOs regions, the linear dependence of the probability and the observables transforms into 
quadratic dependence. 




(27) 



Pee~l-(l-e')sin2| , 



(28) 



7 




AmV(4E) (eV) 

FIG. 1. The dependence of the averaged Vf. survival probabihty inside the sun on AmV4-E for 
different values of e (numbers at the curves). The upper curves describe the corresponding negative 
values of e. Solid lines show results of numerical calculations. Dashed lines correspond to analytical 
formulae with exponential approximation for the solar density profile. 



D. Regeneration effects in the Earth 



For a neutrino arriving at night time, Earth matter effects should be taken into account. 
To gain a qualitative understanding of the Earth effects, we make the crude approximation of 
uniform density, such that the mixing angle in the Earth is constant, Qe-, along the neutrino 
trajectory. In this case, the neutrino propagation has an oscillatory character. We assume 
that the oscillations are averaged out, and therefore 



P2e (vacuum) = sin^ 6*, 



Pae (uniform density) = - [sin^ d + ^\^{2Qe - ^)] • (29) 
It is convenient to introduce a regeneration factor which describes the Earth matter effect: 
/reg = P2e(Hiatter) — P2e(vacuum) = - sin 26'£; sin(26'£; — 26'). (30) 
Denoting by r;^; the parameter r\ in the Earth: 

= vipEyeE), (31) 

where t] is defined in Eq. (13), we find from Eq. (30): 
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Note that /reg is always (for any value tyrr? jE and e) positive, i.e. the matter effect of 
Earth always enhances the survival probability Pee- 
The parameter /reg can be expanded around e = 0: 



/reg 



2(l + r/|) 



1 + 



1 + ^1 



2ve 



) 



(33) 



The regeneration factor has a maximum at r^^; ~ 1 which corresponds to Arn^/E ~ (2 — 
3) X 10~^^ eV. This determines the strong regeneration region in the Am'^/E scale. It is 
situated in the middle of the the solar MSW region. Strong regeneration effects are already 
excluded by the SuperKamiokande result on Day- Night asymmetry [32,33]. Consequently, 
the strong regeneration region separates two (allowed) parts within the MSW region in which 
the regeneration effects are small: 

1. HIGH Am^ region, where r^e 3> 1 and /reg ~ ^/ve', 

2. LOW Am^ region, where tje <^ 1 and /reg ~ ^e/2. 

We will quantify the borders of these regions in Sec. III. In both cases the regeneration 
effect is suppressed by a small parameter and it disappears when moving away from the 
strong regeneration region. 

Let us stress that the SuperKamiokande limit on regeneration effects holds for the energy 
range £■ = 5 — 15 MeV to which this experiment is sensitive. This corresponds to Am^ = 
(2 — 3) X 10~^ eV^. However, strong regeneration effects are not excluded for other energies, 
in particular, for low energy neutrinos. In this case tje can be of the order one and /reg at 
its maximum. 

To first order in e, we obtain from Eq. (33): 



In both the HIGH and LOW regions the dependence of the regeneration factor on e is further 
suppressed by the small parameter imn{2/r)E, ^tje}. 

The probability Pgg at night time, P^, can be written (for the region where the oscillating 
term in Eq. (9) is absent) as follows: 




(HIGH) 
(LOW) 



(34) 




The average daily survival probabihty is given for 775 <S 1 (cos 29s 



l)by 



P ^ ^{Pn + P^) = 1[1 + (1 - 2Pi)(/reg - cos 2^)]. 



(36) 



The Day-Night asymmetry is given by 



Pn — Pd 
P 



2/reg 



(37) 



l/(l-2Pi) -cos2^ + /reg' 
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Let us now consider the dependence of An-d and P on e in the HIGH and in the LOW 
regions keeping the lowest order terms in e and in mm{riE, I/ije}- 

1. In the HIGH region, the adiabaticity condition is satisfied and we can safely put 
Pc — 0. Consequently, Pi = cos^ 9s- Then, Eq. (36) simphfies: 



P=-[l-cos2^5(/reg-cos2^)]. 



(38) 



With near-maximal mixing, we obtain: 



^ 1 / _ e- l/2rjE 



1 + 



(39) 



For ry| -C 1 (small Am^) we find: 



2 V 277E r)i;J 



(40) 



where the last term is a small correction which comes from the regeneration factor. For the 
Day-Night asymmetry we obtain: 



1 



'II 



Vs 



1+^5 



(41) 



Notice that the e-dependent effect changes sign between large and small rjs- For 77I ^ 1 we 
get: 



An-d — (1 + e) 
Ve 



(42) 



The asymmetry increases linearly with e for small enough rjg. 

2. In the LOW region we can take rjs — and, consequently. Pi = Pc- Then, Eq. (36) 
simplifies: 



1 r 



l + (l-2Pe)(/reg-e) 



(43) 



In LOW region we have 7^1 (see Eq. (23)) and, therefore, Pc <C 1 (see Eq. (22)). In any 
case, the Pc-dependent term in Eq. (43) is suppressed by a small factor and can be neglected. 
We get: 



P^l{l-e + 



r]E_ 
2 

For the Day- Night asymmetry we obtain taking Pc <S 1: 

^ _ , 2/reg 



IN-D 



l + 2Pc-e + /reg' 



(44) 



(45) 



and in the limit t^e <^ 1: 
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An-d ^Ve{i- 2e-^/' - Y + ■ (46) 

A few comments are in order. The Day-Night asymmetry is suppressed by tje- The e- 
dependent effect is one of several small corrections to the leading result. For large parts of the 
LOW region it is the leading correction, but for large Am'^/E, the subleading regeneration 
effect could be comparable while for small Arri^/E, the non-adiabatic correction could give 
the main correction. 

As we mentioned before, a strong regeneration effect with ~ 1 and /rcg --^ 1 is not 
excluded for low energy neutrinos. In particular in the LOW region strong regeneration can 
show up for the beryllium and pp-neutrino components of the spectrum. In this case, the 
approximation tje <^ 1 does not work and one should use the complete expression for the 
regeneration factor and the asymmetry. 

As seen from Eqs. (40), (44), (42) and (46), the dominant dependence on e of both P 
and ^N-D arises from the dependence of the solar survival probability on e. The dependence 
which follows from the regeneration factor is further supressed by min{r;£;, I/tje}- (For A^.d 
the e dependence follows from the dependence on P in the denominator.) 

From Eqs. (42) and (46) we find that the Day-Night asymmetry strongly depends on 
Am^ due to the rjE or I/tje factors. In the HIGH region An-d oc Am^ while in LOW region 
An-d oc 1/Am^. The dependence of the asymmetry on e is much weaker: In both regions 
^N-D oc (1 + e). Consequently, the measurements of ^n-d are very sensitive to Am^ while 
the sensitivity to e is substantially lower. 

Let us comment on the range of validity of the approximate treatment of the Earth 
effects [25]. In the HIGH region, for E ~ MeV and Am^ ^ lO^^eV^, the oscillation length 
li, is much shorter than the size of the Earth and neutrinos undergo many oscillations inside 
the Earth. The constant electron number density approximation gives a good description 
of An_d which involves integration over the zenith angle. In the LOW region, for Am^ <S 
10~^eV^, the oscillation length is approximately equal to the refraction length Iq and the 
latter is comparable to the size of the Earth (independently of Am^). The details of the Ne 
profile do not play an important role, and the effect is determined by the average density 
along the neutrino trajectory. 

Regeneration effect leads also to seasonal variations of signals [25-27]. These variations, 
however, are less sensitive to the oscillation parameters. 



E. Distortion of the energy spectrum 

The distortion of the solar neutrino energy spectrum can be characterized by the distor- 
tion parameter defined as: 

^dP 

'^^PdE- ^''^ 

Averaged over the appropriate energy interval this parameter is proportional to the shift of 
the first moment of the spectrum or to the slope parameter used in the literature. 

To understand the distortion of the spectrum (energy dependence of the averaged prob- 
abihty), we remind the reader that r] oc 1/E and use Eq. (38) for the HIGH and QVOl 
regions and Eq. (43) for the LOW region. 
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1. The HIGH and QVOl regions: 

(i) For large Am^, the effect of the adiabatic edge of the suppression pit which is encoded 
in the rjs dependence is important: 

The distortion parameter is proportional to e. The slope (shift of the first moment) is 
positive (negative) for e < (e > 0). The distortion decreases rapidly with rjs- 

(ii) For small Am^, the earth regeneration which is related to the r]E dependence is 
important: 

1 +iiV (49) 




277b "He, 

The distortion decreases with increase of 77^ or Am^. The sensitivity to e is much weaker than 
in the previous case and it follows mainly from the dependence on the average probability 

in the denominator. 
2. The LOW region: 

(i) For large Am^, the regeneration effect is important: 

«.«^f(l + 6). (50) 

The slope increases with rjE- The dependence on e is weak. 

(ii) For small Am^, the effect of the nonadiabatic edge of the solar suppression pit gives 
the dominant effect: 

rIP 

= 2eE^ ^ e^e-^l\ (51) 

OjEj 

The distortion is proportional to e. The slope is positive (negative) for e > (e < 0). The 
effect is suppressed for relatively weak violation of the adiabaticity. 

As we have mentioned previously, for non-averaged vacuum oscillations s^^ oc (1 — e^), 
that is, the dependence of distortion on is very weak. 



F. Summary 

Let us summarize the results of our analytical studies. We have found simple analytical 
expressions for various observables (rates of events given by the survival probability, Day- 
Night asymmetry, distortion of the spectrum etc.) in terms of Am^ and the deviation 
from maximal mixing e. These approximate analytical expressions reproduce correctly the 
functional dependence of the observables on Am^ and e and allow us to understand all the 
features of the exact numerical calculations. 

The effects and the dependence of observables on e change with Am^. Accordingly, 
we define several regions of Am^ which have different physical pictures. As Am^ decreases 
from its upper (CHOOZ) bound, we have the following regions (we quantify borders of these 
regions in the next section): 



12 



• quasi- vacuum oscillation region with large Am^ (QVOl), 



• MSW region with high Am^ (HIGH), 

• MSW region with low Am^ (LOW), 

• quasi- vacuum oscillation region with small Am^ (QVOs), 

• region of non-averaged vacuum oscillations (VO). The high Am^ part of the VO region 
will be called VACl- 

As concerns the e dependence of observables in these regions, we find two main conclu- 
sions: 

1. Maximal mixing, e = 0, is not a special point as far as the phenomenology is concerned 
(in contrast with theory). No divergencies or discontinuities appear in the dependence of 
observables on e. The dependence on e is smooth and, in many cases, very weak. To mention 
a few examples: The day time survival probability is 1/2 at e = 0. Earth regeneration ef- 
fects, however, enhance the survival probability. In the MSW region the slope of the energy 
spectrum distortion is proportional to e in certain regions of Am^, and consequently changes 
sign at e = 0. The Day- Night asymmetry is proportional to (1 -|- e). 

2. The character of the e-dependence of observables is different in the vacuum oscilla- 
tions and MSW regions. In the regions of vacuum oscillation all the effects depend on e 
quadratically. More precisely, the e-dependent factors are of two types: 

1 + (averaged oscillations), 

1 — (non — averaged oscillations). (52) 

The dependence is symmetric with respect to interchanging e <-> — e. In the MSW regions 
(both HIGH and LOW), observables depend on e linearly. Obviously the dependence is non- 
symmetric with respect to e ^ — e. In the region of strong adiabatic conversion (bottom 
of the pit) we get for the survival probaility, the Day-Night asymmetry and the distortion 
parameter: 

Pd oc (1 - e), An-d, s„ (X (1 + e). (53) 

At the edges of the MSW region (edges of the suppression pit) we find 

Su oc e. (54) 

The quasi- vacuum oscillation regions, QVOl and QVOs, are transition regions between the 
MSW and vacuum oscillation regions, where the linear dependence of observables transform 
into quadratic dependence. 

Thus, for e ^ 1 the sensitivity of experiments to deviation from maximal mixing is much 
higher in the MSW regions. It will be difficult to measure e near maximal mixing if Am^ 
turns out to be in the VO or QVO regions. 
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III. PRESENT STATUS OF MAXIMAL MIXING 



In this section we describe the present status of maximal as well as near-maximal mixing. 
We use for this purpose the latest available results on solar neutrinos from Homestake [28] , 
SAGE [29], GALLEX+GNO [30,31], and the 1117 days of data sample of Super-Kamiokande 
[32,33] experiments. 

We calculate the acceptability of maximal and near-maximal mixing as a function of 
Am^ in the whole allowed range, that is, below the upper bound Am^ < 10~^ eV^ from 
the reactor experiments [16,34]. The goal of this study is to find excluded regions of the 
oscillation parameters Am^ — e, as well as the regions of these parameters which are allowed 
and most plausible. Wc quantify our statements in terms of confidence level at which a given 
region is accepted (probability of realization) or excluded. Wc study the dependence of our 
conclusions on uncertainties in the solar neutrino fiuxes (the SSM uncertainties) as well as 
on the procedure employed in the analysis. Some model- and the procedure-independent 
statements are formulated. The extent to which the results of this Section hold in a three 
generation framework with \Ue3\ 7^ is discussed in Sec. VI. 

A. Global fit: allowed and forbidden regions 

The results of a global fit to all existing experimental data on solar neutrinos are shown 
in Figs. 2—5. The analysis includes rates in Chlorine [28], GaUium [29-31], and Super- 
Kamiokande [32,33] experiments, as well as the zenith angle dependence [32,33] and the 
shape of the recoil electron spectrum [32.33] in Super-Kamiokande. 

In Fig. 2 wc plot the contours of constant confidence level (iso-contours) in the Am^ — 
e plane. Points inside a given contour are accepted at a lower confidence level than on 
the contour itself. In the "global" analysis we combine the information on the Day-Night 
variation of the event rates and the recoil energy spectrum at Super-Kamiokande by using 
their independently measured spectra during the day and during the night. With this the 
total number of independent experimental inputs in the global analysis is 38 which includes 
3 rates, and 35 data points for the Super-Kamiokande day and night recoil energy spectra 
(2 X 18 bins minus 1 overall normalization). We do not include in the analysis the new 
lower energy bin as its systematic uncertainty is still under study by the Super-Kamiokande 
Collaboration [33]. We use the solar neutrino fiuxes from the Standard Solar Model (SSM) 
of Ref. [35] (BP98). The contours have been defined by the shift in x^, Ax^, with respect 
to the global minimum in the whole plane of the oscillation parameters. The minimum lies 
in the LMA solution region: 

= 33.4 for 36 d.o.f. (55) 

which corresponds to a probability of 59%. For details on the statistical analysis we refer to 
Ref. [4]. 

In Fig. 3 we show the dependence of A^^ on Am^ for three different values of e 
(— 0.3, 0, +0.3). This figure corresponds to three Am^-profiles (cuts) of the confidence level 
from Fig. 2. 

What is the impact of individual experimental results on the global fit? Panels (a) — (f) of 
Fig. 4 show the contours of constant for each individual observable i. As mentioned above 
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FIG. 2. Contours of constant confidence level (iso-contours) in the Am^ — e plane. The shaded 
regions are accepted at the given or lower confidence level by the global fit of all available solar 
neutrino data. The white area is excluded at 99.99% confidence level. 

the total in panel (f ) is obtained by combining the of the individual rates (including the 
correlation of their theoretical errors) with the corresponding for the Super-Kamiokande 
night and day recoil energy spectra. 

Our calculations allow us to define regions of oscillation parameters that are excluded 
and other that are allowed. It follows from Fig. 2 that there are two main regions of Am^ 
which are excluded at a very high (more than 99.99%) confidence level: 

1. The regeneration region: for maximal mixing we find 

Am^ = (0.6 - 8) X 10"^ eV^ (e = 0). (56) 

The excluded region increases with e for positive values of e. In particular, for e = 0.3, the 
excluded range becomes Am^ ~ (0.4 - 10) x 10"^ eV^ The excluded region also expands 
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Am' (eV') 

FIG. 3. The dependence of Ax^ from the global analysis (with the boron neutrino flux of the 
SSM) on Am^ for e = —0.3 (dashed line), (solid line) and +0.3 (dash-dotted line). Dotted 
horizontal line marks 99% CL. Below this line the corresponding oscillation parameters are accepted 
at a confidence level lower than 99%. 



for negative values of e at e <— 0.25. 

In the regeneration region the solar neutrino observables are strongly modified by the 
z/g-regeneration in Earth. As we have discussed in Sec. II, the regeneration is always pos- 
itive, thus leading to an increase in the z/g-flux. Correspondingly, the counting rates in all 
experiments increase. There are two main physical effects in the regeneration region that 
are inconsistent with observations and therefore lead to the exclusion: 

(i) A large Day-Night asymmetry, An_d ^0.2 (see Figs. 4(d) and 6), is in contradiction 
to the Super-Kamiokande result and plays a dominant role for positive e. 

(ii) A large Ar-production rate, Qat > 4 SNU (see Figs. 4(b) and 6), is in contradiction 
to the Homestake result and leads to the increase of the excluded region for negative e. 

2. The vacuum oscillation region: 

Am^ = (0.14 - 4) X 10"^° eV^ (e = 0). (57) 

The size of this region depends very weakly on e in the interval —0.3 < e < +0.3. The 
exclusion follows from the interplay between the total rates and the shape of the recoil 
electron energy spectrum. Notice that the rates only (see Fig. 5(a)) can be reproduced 
rather well in some parts of the region in Eq. (57). However, in these regions the distortion 
of the recoil electron spectrum is in contradiction with the Super-Kamiokande results (see 
Fig. 4(e)). More specifically, for Am^ ~ 10~^° eV^, a negative slope (or shift of the first 
moment) in the "reduced" spectrum is expected. The reduced spectrum is defined as the 
ratio 

R{E) ^ N,{E)/N,{Ef^'^, (58) 
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1 nr^^ r — 1 — I I I I I I I I I I I I I 

0,3 0,2 0.1 0.0 -0,1 -0,2 -0,3 0,2 0,1 0,0 -0.1 -0,2 -0,3 

FIG. 4. Contours of constant in the Am^ — e plane for individual experimental results. The 
various panels show the contributions of different results to the total 'x^: (a) total rates in SAGE 
and GALLEX, (b) total rate in Homestake, and the Super-Kamiokande measurements of (c) the 
total rate, (d) zenith angle dependence, and (e) recoil electron energy spectrum. In (f) we show 
iso-contours of the total \^ . 

where Ni{E) {Ni{E)^^^) is the number of events in a given electron energy, E, bin i with 
(without) oscillations. In the same way we define any "reduced" observable as the ratio of 
its value with respect to the expected one in the SSM in the absence of oscillations. 

3. There is another region forbidden at 99.99 % CL extending from Am^ = (0.3 — 2) x 
10~^ eV^ for e = 0.3 to Am? = 6 x 10~^ eV^ and e = 0.22. As seen in Fig. 3 this region is 
forbbiden with a lower CL. The exclusion follows mainly from the effect of the rates as seen 
in Fig. 5(a) being mainly driven by the bad fit to the Gallium rate (see Fig. 7(a)). 
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We distinguish five regions of tlie oscillation parameters where maximal mixing is allowed 
at a confidence level that is lower than 99.9%: 

1) Quasi vacuum oscillation region with large Am^ (QVOl): 

where E is the average detected energy for a given experiment. The upper bound comes 
from reactor experiments [16,34]. Here the flavour conversion is mainly due to averaged 
vacuum oscillations with only small matter corrections inside the Sun and the Earth. 

2) MSW region with high Am^ (HIGH): 

This region corresponds to the maximal and near-maximal mixing part of the LMA solution. 
It is restricted from below by strong Earth regeneration effects (large Day-Night asymmetry 
and large Ar-production rate). Maximal mixing is acceptable at confidence level larger than 
99.1%. As follows from Fig. 3, the dependence of Ax^ on Am^ is rather weak. The global fit 
becomes substantially better with increase of e (shift to positive values): e > 0.25 is accepted 
at 90% CL. For e = 0.3, the 90% CL allowed region expands to Am^ = (2 - 10) x 10"^ eV^. 
In contrast, the goodness of the fit decreases when we shift to negative values of e. 

3) MSW region with low Am^ (LOW): 

Am^ ~ (0.1 - 3) X 10"^ eVl (61) 
Here maximal mixing is acceptable at ~ 99% CL in the interval 

Am^ - (0.1 - 2) X 10"^ eV^ (CL < 99%). (62) 

For positive e the fit improves while for negative e it worsens. In particular, e = 0.2 is 

accepted at ~ 90% CL for Am^ ~ (0.8 — 1.5) x 10~^ eV^. The local minimum occurs at 
Am^ = 1.0 X 10"'' cV^ and e = 0.21. With increase of e the accepted region of Am^ shifts 
to larger values. At e = 0.3 we obtain Am^ = (0.6 — 2) x 10"'' eV^. Conversely negative 
values of e are disfavored. 

4) Quasi vacuum oscillation region with small Am^ (QVOs): 

Am^ - (0.1 - 1) X 10"^ eVl (63) 

In this region the flavour conversion is due to (mainly non-averaged) vacuum oscillations 
with small matter effects. Maximal mixing is acceptable at ~ 99% CL in the interval 

Am^ ~ (0.5 - 1) X 10-^ eV^ (CL ~ 99%). (64) 

The flt becomes worse with increase of |e|, but while for e > the QVOs region is essentially 
excluded, for e < we still have a reasonably good flt. 

5) Vacuum oscillation region with relatively large Am^ (VACl): 

Am^ ~ (0.4 - 1) X 10"^ eV^ (65) 
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Maximal mixing is accepted at a confidence level better than 99% only in a very small 
interval centered at: 

Am^ ~ 6.6 X 10-^° eV^ (CL < 99%). (66) 

In this interval, the goodness of the fit depends on e very weakly. 

Summarizing, maximal (or near-maximal) mixing is allowed at 99% or slightly lower CL 
in several small intervals of Am^ in the QVOl, HIGH, LOW, QVOs and VACl solution 
domains. The values e = 0.05, 0.1, 0.2 are allowed at 99, 95 and 90% CL, respectively. At 
4(7, practicahy the whole HIGH, LOW, QVO and VACl ranges are allowed. 

Let us point out the role of individual experimental results in constraining maximal 
mixing (see Fig. 4). The rates in the Gallium and Super-Kamiokande experiments can be 
well accounted for at maximal (or near-maximal) mixing, although the Super-Kamiokande 
measurement slightly disfavors a negative e. The zenith angle distribution measured by 
Super-Kamiokande gives some preference to the HIGH region and excludes the strong re- 
generation region. In contrast, the Super-Kamiokande result on the recoil electron energy 
spectrum gives some preference to the VACl, HIGH and LOW regions and excludes the 
range Am^ = (0.3 — 4) x 10~^° eV^. Both the zenith angle distribution and the shape of 
spectrum have weak dependence on e. In contrast, total rates are sensitive to e, especially 
in the HIGH and LOW regions. 

B. Dependence of the results on features of the analysis 

Let us study the dependence of the allowed and excluded regions in the Am^ — e plane 
on features of our analysis. 

Total rates versus spectrum and zenith angle distribution: The total rates give the most 
stable, reliable, and statistically significant information. Wc have carried out a fit to the 
three rates only. In Fig. 5(a) we show the dependence of the shift of for this analysis with 
respect to the absolute minimum in the whole plane of oscillation parameters. The absolute 
minimum for the analysis of the three rates, — 0-76 (for one d.o.f) is achieved in the SMA 
region. Comparing Fig. 5(a) with Fig. 3, we learn that if we exclude the information from 
the recoil electron energy spectrum and the Day-Night variation of the event rates from our 
analysis, then the allowed and forbidden regions are substantially modified. In particular, 
we note in Fig. 5(a) the following three features: 

(i) The goodness of the fit at maximal mixing from the analysis of the three rates only is 
worse in whole MSW region. At 99% CL only small interval in the LOW region is allowed. 
The fit improves however with increase of e. This means that it is the data on the spectrum 
and the zenith angle distribution which favor maximal mixing. 

(ii) Allowed regions appear in the VAC solution range. Wc learn that the data on the 
spectrum and the zenith angle distribution exclude (otherwise) allowed VAC regions. 

(iii) The regeneration region is still strongly disfavored by the high Ar-production rate. 

The Homestake result: Consider the impact of the Ar-production rate on our results. 
In Fig. 4(b) wc show the fit to only this rate. From the figure we see that the Homestake 
result strongly disfavors maximal mixing for all Am^ above 10^^° eV^, that is, in all the 
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globally allowed regions. In Fig. 5(b) we show the result of a global fit to the data without 
the Homestake result. Clearly, the acceptability of maximal mixing improves for all Am^ > 
10~^° eV^ with the best fit points being in the LOW region, Am^ = 10~^ eV^. For e = 0.3 
the best fit is in VACl region. We would like to emphasize the following points about a fit 
without the Homestake result: 

(i) In the HIGH region, maximal mixing is accepted already at 1.7 cr with very little 
dependence on e. 

(ii) In the whole LOW region, maximal mixing gives a very good fit. 

(iii) In the whole QVOl region, maximal mixing gives a very good fit, but positive values 
of e are still disfavored. 

(iv) Regions of strong regeneration and VAC (small Am^) solutions are excluded by the 
Super-Kamiokande data on the spectrum and the Day-Night asymmetry. 

Solar neutrino flux uncertainties: Of all the relevant solar neutrino fiuxes, the boron neu- 
trino fiux suffers from the largest uncertainty, leading to systematic errors in the predicted 
detection rate that cannot be estimated reliably at present. One way to avoid this problem 
is to determine the boron neutrino flux experimentally, using the total rate measured in the 
Super-Kamiokande experiment. Similar results are obtained if the boron neutrino flux is 
treated as a free parameter in the analysis (in this case the Super-Kamiokande rate, being 
the most precise and sensitive to the boron neutrino fiux, will fix this fiux). In Fig. 5(c) we 
show the results of a global fit with the boron neutrino fiux treated as a free parameter. We 
plot the dependence of Ax^, the shift with respect to the absolute minimum, on Am^ for 
three different values of e. The shape of the curves is very similar to that in Fig. 3. Also the 
allowed and excluded regions at a given CL practically coincide with those in Fig. 3. (One 
should take into account that now we have one additional free parameter and therefore a 
99% CL corresponds to Ax^ = 11.36). The reason of this similarity hes on the fact that 
both the spectrum and the Day-Night variation are flux independent. Furthermore in both 
cases the boron neutrino flux is flxcd by the Super-Kamiokande result. 

A comparison of Fig. 3 and Fig. 5(c) shows that the results of our analysis are stable 
with respect to the way that the uncertainty in the boron neutrino flux is treated. 

IV. MAXIMAL MIXING AND PREDICTIONS FOR INDIVIDUAL 

EXPERIMENTS 

In this section we consider the predictions of (near-) maximal mixing for various observ- 
ablcs. This will clarify the sensitivity of individual experiments to the neutrino oscillation 
parameters in the relevant range. The extent to which the results of this Section hold in a 
three generation framework with |?7e3| 7^ is discussed in Sec. VI. 

A. Total rates 

In Fig. 7 we plot the values of the expected event rates: Germanium production rate. 
Argon production rate, and the rate of events at Super-Kamiokande as functions of Am^ 
for 3 values of e: +0.3, 0, —0.3. The rates are normalized to the no oscillation expectation, 
i?j = Qi/Qf^^. For each value of e we plot three curves: the central curves give the expected 
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rates using central vahics of the BP98 fluxes and the Tipper and lower lines represent the 
theoretical uncertainty (without the error for the interaction cross sections) from varying the 
nine parameters in the SSM within ilcr. The horizontal lines give the experimental values 
within ±1(7 experimental errors. The vertical lines in the range Am^ = (3 — 10) x 10~^ 
eV^ give the expectation from the SMA 99% CL region (again, including the theoretical 
uncertainties in each point). 

Let us first discuss the dependence of the rates on Am^. The rates are proportional to 
the survival probability. Therefore, the main features of Fig. 7 reflect the dependence of 
the survival probability on iSm^ jE (see Fig. 1 and Eqs. (26), (27) and (28)). For maximal 
mixing the survival probability in the Sun as a function of iS.vr? jE is constant: Po — 1/2 
(see Eq. (19)). The probability Pee is enhanced by the Earth regeneration effect for !\vn^ jE 
in the range (10~^^ — 10~^^) eV (sec Eq. (36)). For e 7^ 0, the effects of the adiabatic edge 
situated at Am^/E" = (10~^^ — 10"^°) eV (Eq. (40)), and of the non-adiabatic edge situated 
at ^m^jE = (10-^^ - 3 X 10"^^) eV (Eq. (43)), become important. For ^m^jE ^ lO'^^ eV, 
an oscillatory behaviour appears due to non-averaged vacuum oscillations between the Sun 
and the Earth (Eq. (28)). With certain modifications all these features can be seen in Fig. 7. 
The simplest dependence is for the Super-Kamiokande rate since only one (boron) neutrino 
flux contributes (Fig. 7(c)). The Ar-production rate has an additional fine structure due 
to contributions from additional fluxes, and in particular the beryllium neutrino flux. As 
can be seen from Fig. 7(b), an additional enhancement appears in the regeneration region 
at Am^ = 3 X 10~^ eV^ and the probability as a function of Am^ becomes asymmetric 
in the regeneration region. In the VAC and QVO^ regions the boron neutrino "wave" is 
modulated by the beryllium wave with a smaller amplitude. For the Ge-production rate, 
all the features of the curves are shifted (with respect to the Super-Kamiokande curves) 
to smaller values of Am^ by a factor ~ 30. This feature is due the fact that the main 
contribution to Qq(, comes from the pp-neutrino flux with an average detected energy of 0.3 
McV, about 30 times smaller than the average energy of the boron neutrino flux. 

Let us consider now the e-dependence of the rates. We distinguish here between three 
different regions of Am^: 

1) The QVOl region with Am^ ^ 3 x 10~^ eV^. Here we have essentially averaged 
vacuum oscillations, so that the survival probability is given by Eq. (26). 

2) The QVOs and VAC regions with Am^ < 10~^ eV^. Here we have essentially non- 
averaged or partially averaged vacuum oscillations which take the form of Eq. (28). In 
principle matter effects strongly suppress the oscillations inside the Sun and the Earth. 
However, the modification of the observables is small, since the size of the Sun (and the 
Earth) is much smaller than the oscillation length in vacuum. 

3) The matter conversion (MSW) region is between these two QVO regions. The oscilla- 
tion effect is strongly suppressed. As explained in Sec. II, Earth matter effect is important 
but insensitive to e. The expression of P in this region is given in Eqs. (40) and (43). 

From Eqs. (26), (27), and (28) we can find straightforwardly the dependence of the 
observables on deviations from maximal mixing. For the VO, QVOs ^-nd QVOl regions, the 
following points are in order: 

1. The observables depend very weakly on e. Corrections to maximal mixing are of C(e^). 

2. The dependence is symmetric under the exchange e ^ — e. Minimum of the survival 
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probability, and consequently, minima of rates, are at e = 0, that is, at maximal 
mixing. 

For the MSW regions, the following points are in order: 

1. The survival probability and the rates depend linearly on e. Corrections to the maximal 
mixing case are consequently larger. 

2. The survival probability and the rates decrease with increase of e. 

There are two transition regions between VO and pure matter conversion. In these 
regions, the symmetric dependence of the observables transforms into a linear dependence 
and the sensitivity to deviation from maximal mixing increases. The ambiguity e — e 
disappears. 

Thus, in the region of pure matter conversion (inside the Sun) the sensitivity of measure- 
ments of rates to a deviation from maximal mixing is maximal. It is in this region that the 
possibility of maximal mixing can be tested with the highest accuracy. The corresponding 
Am^ range depends on the neutrino energy. For the highest energies of the solar neutrinos 
{E ~ 10 MeV) we get the range of maximal sensitivity: Am^ = 10~^ — 10~^ eV^, while for 
the lowest energies (£■ ~ 0.3 MeV) the range is Am^ = 3 X 10-9 - 3 X 10"^ eV^ Studying 
effects by the experiments with different energy thresholds we can get high sensitivity to 
deviations from maximal mixing in whole range of Am^ excluding VO. 

In the next two subsections we consider the dependence of specific rates on the oscillation 
parameters and evaluate the sensitivity of their present measurements to deviations from 
maximal mixing. 

B. Ar— production rate. 

Let us now study the implications of the results presented in Fig. 7(b) for the Ar- 
production rate and in Fig. 7(c) for the Supcr-Kamiokande rate in the various Am^ regions. 
As can be seen in Fig. 7(b), the Ar-production rate for maximal mixing in all favorable 
regions (HIGH, LOW, VACl) lies in the range Rci = 0.50±0.08 {Qat = 3.9±0.6 SNU), which 
is 2a above the Homestake result. The predictions for e = +0.3 (—0.3) are below (above) the 
values for maximal mixing by just about la. As mentioned above, the highest sensitivity 
of Qat to e is achieved already in the HIGH region, where Qat depends on e linearly. 
For the Super-Kamiokande rate, we get that for maximal mixing in all favorable regions 
RsK — 0.58 ± 0.12 which is la above the Super-Kamiokande result. Moreover comparing 
Figs. 7(b) and 7(c) we see that the Ar-production rate and the Super-Kamiokande rate 
are strongly correlated as they are both dominated by the contribution from boron flux 
neutrinos. 

Next, in order to reduce the SSM uncertainty we normalize the boron neutrino flux to 
the Super-Kamiokande rate. More precisely, for each pair of the of oscillation parameters 
(Am^, e) we flnd the boron neutrino flux which reproduces the Super-Kamiokande event 
rate. All other fluxes and their uncertainties are taken according to the BP98. Using this 
procedure we calculate Qat (and in the next subsection also Qoe)- The results of this 
calculation are shown in Fig. 8. 
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Prom Fig. 8 we see that after boron flux normalization procedure the dependence of Qat 
on e is relatively weak since the ratio between the boron neutrino flux contribution to Q^r 
and the contribution from charged current interactions to the Super-Kamiokande rate is 
independent of the survival probability and therefore of e. The e dependence of Qat comes 
directly from the suppression of the beryllium and other (CNO and pep) neutrino fluxes, and 
indirectly from the contribution of neutral current interactions to the Super-Kamiokande 
rate. Expressing the boron neutrino flux Jb via the SK rate: 

P{1 — r) + r 

we obtain: 

Here r is the ratio of the f ^ — e and — & cross-sections, Q^^ and Q^r ^-^e the SSM 
contributions to the Ar-production rate from the boron neutrino flux and from all other 
low energy fluxes, respectively, and P and P' are the effective survival probabilities for the 
boron neutrino and low energy neutrino fluxes, respectively. The bands in Fig. 8 reflect la 
errors due to the uncertainties in all, but the boron, neutrino fluxes. 

As expected, in the QVO regions, QVOl and QVOs, the e-dependence is symmetric 
around e = 0, Qat — QAr(e^), while in the pure matter conversion regions, HIGH and LOW, 
Qat depends on e linearly. In the linear regime, the change in the Ar-production rate is 
^Qat- ~ (0.7 — 0.8) SNU for —0.3 < e < +0.3. This is about 3(T for the present experimental 
error. In the quadratic regime the change is substantially smaller: ^Qat ^ (0.1 — 0.2) SNU. 
Clearly, the present sensitivity is not enough to draw deflnite conclusions. Moreover, even 
after normalization of the boron neutrino flux to the Super-Kamiokande rate, the predicted 
rate is higher than the Homestake result for all globally allowed values of Am^. The only 
statement that one can make is that the Homestake result favors a significant deviation 
from maximal mixing in the MSW region. Therefore checks of the Homestake result and 
improved accuracy in measurements of Qat by a factor of two or higher would have important 
implications for maximal and near-maximal mixing. 



C. Ge— production rate 

Let us now study the implications of the results presented in Fig. 7(a) where we show 
the Ge-production rate as a function of Am^ for different values of the e parameter and 
in and Fig. 9 where we plot Qoe as a function of e in the various Am^ regions. While 
the results plotted in Fig. 7(a) are obtained within the SSM, the predicted rates shown 
in Fig. 9 have been obtained after normalization of the boron neutrino flux to the Super- 
Kamiokande measured rate as described previously for Fig. 8. In this case, unlike in the case 
of Ar-production rate, the results are very slightly modifled by the flux normalization 
since the Ge-production rate is dominated by the contribution from the pp neutrino flux 
and the corresponding theoretical uncertainties are smaller. 

In the QVOl region, the rate is determined by averaged vacuum oscillations and therefore 
the expected rate is symmetric around maximal mixing. According to Eq. (28), P{e — 
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±0.3) = 0.545 which corresponds to Qce = 70.3 ± 3.5 SNU (with la theoretical error). For 
exactly maximal mixing we have P(e = 0) = 0.5 and the rate is minimal, Qce = 64.5 ± 3.5 
SNU. As seen in Fig. 7(a) the bands for the different e-values almost overlap and all the 
predictions are in agreement with the present experimental result. In other words, with 
present experimental error bars it is impossible to measure deviations from maximal mixing. 

In the part of the HIGH region with small Am^ the effects for pp neutrinos occur in the 
transition between linear and quadratic e regimes. Consequently, the sensitivity of Qce to 
e is still below the present experimental accuracy. We find from Fig. 7(a): Qcei^ = 0) ~ 
g^g(e = +0.3) 60 - 68 SNU, while Qcd^ = -0.3) 70 ± 4 SNU, approximately la 
(experimental) higher. From Fig. 9(a) we see that after the normalization the variation 
of the Qce in e range (—0.3, +0.3) is AQ^e ~ (10 — 12) SNU, which is about 2a (present 
experimental error). 

In the LOW region the sensitivity to e is maximal: the pp neutrinos undergo pure matter 
conversion and the rate depends on e linearly. We get, for Am^ = 10~^ eV^: Qcei^ — 
-0.3) = 93 ± 4 SNU, Qoeie = 0) = 77.5 ± 3.5 SNU and Qceie = 0.3) = 61 ± 3 SNU. The 
difference between the predicted values of Qcei^ = —0.3) and Qcei^ = +0.3) is more than 
2a (experimental). From Fig. 9(b) we see that after the normalization the variation of 
the Qce in £ range (—0.3, +0.3) is AQce ~ (33 — 35) SNU, which is more than 5a (present 
experimental error). In this case one gets certain restrictions on e, although the confidence 
level is low. For example, for Am^ = 10"^ eV^ the combined SAGE and GALLEX+GNO 
result gives the la range —0.12 < e < +0.2. Therefore, further decrease of the experimental 
error bars by a factor of two, from the present 5 SNU to 3 - 4 SNU, could have important 
implications for the mixing provided that Am^ will be fixed by some other independent 
measurement. Notice that in the LOW region one expects maximal regeneration effect for 
the pp-neutrinos which can be detected as, e.g., seasonal variation of the Ge-production 
rate [27]. 

The situation is similar in the QVOs region down to Am^ = 5 x 10~^ eV^: the Ge- 
production rate depends on e linearly and AQce ~ 30 SNU (sec Fig. 9(c)). 

In the VO region, deviations from the maximal mixing result are determined by and 
the variations (for a given Am^) are small as shown in Fig. 9(d): AQoe ~ 5 — 7 SNU. Thus 
the sensitivity of present data is still low and practically the whole interval e e (—0.3, +0.3) 
is allowed at (2 — 3)o" level. 

In consequence serious implications for maximal mixing require further decrease of the 
experimental error bars down to 3-4 SNU. 



D. The Day— Night asymmetry in electron scattering events 



In Fig. 6 we show contours of constant Day-Night asymmetry of the u — e scattering 
events in the Am^ — e plane. The Super-Kamiokande 1117 days result 

N - D 

^N-D = 2 ^^^ ~ 0.034 ± 0.026 (69) 

excludes at the 3a level the region An-d ^0.11 which corresponds, at maximal mixing, to 

4 X 10"^ ^ Am^ ^ 1.3 X 10"^ eVl (70) 
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The exclusion interval increases slightly with e. The preferable regions of Am^ for e = are 

Am^ = (2.5 - 10) X 10-^ eV^ Am^ = (0.6 - 2.2) x 10"^ eVl (71) 

We emphasize that these results are SSM independent and have no ambiguities related to 
the analysis of the data. 

E. Zenith angle distribution of electron scattering events 

In Fig. 10 we show the zenith angle distribution of events in Super-Kamiokande for 
maximal and near-maximal mixing and for various vahics of Am^ from the allowed regions. 
Significant enhancement of the night rate is expected in the HIGH and LOW regions. 

In the HIGH region the distribution of events during the night is rather flat and the 
dependence on e is weak, so it will be difficult to use the shape to measure e. The dependence 
on Am^ is also weak. Maximal signal is expected in the N3 or/and N5 (core) bins. 

In the LOW region the oscillations occur in the matter dominated regime (see Sec. II) 
when the oscillation length practically coincides with the refraction length, 1^ ^ Iq- For those 
trajectories crossing the mantle only (N1-N4), the latter can be approximately determined 
by the average density along the trajectory. Maximal effect (which corresponds to the 
oscillation phase = tt) is reahzed for the trajectories with cos9z — 0.3, i.e. in the N2 bin 
(see Fig. 10(b)). The phase 27r is collected at cos^^^ ~ 0.5 which corresponds to a minimum 
of the signal. Notice also that the zenith angle distribution depends on e very weakly. For 
a given trajectory the amplitude of the oscillation is determined by the mixing angle in the 
Earth matter 

sin^ 20E 5 T (72) 

where r]E the parameter r/ in the Earth defined in Eq. (31). Thus sin^2^^ is an increasing 
function of e (see also Eq. (17)). As a consequence the number of events in maxima increases 
with e as seen in Fig. 10(b). Present data do not show any enhancement in the N2 bin. 

We conclude that precise measurements of the zenith angle distribution would allow the 
determination of Am^ and probably resolve the HIGH/LOW ambiguity but are unlikely to 
play a significant role in the determination of e. 



F. The recoil electron energy spectrum 

In Fig. 11 we show the expected recoil electron spectrum for maximal mixing with various 
values of Am^ in the regions allowed by the global fit. In the HIGH and LOW regions the 
"reduced" spectra, R{E), are flat. Strong distortion is expected in the VACl region. Thus 
further improvement on the measurement of the recoil electron spectrum can discriminate 
between the MSW and VACl regions. 

For large enough Am^ and e 7^ a distortion in the spectrum appears due to the effect 
of the adiabatic edge (Fig. 12(a)). This feature can be traced from the dependence of the 
survival probability on Aw? / E (see Fig. 1). For a positive e conversion inside the sun leads 
to negative slope of the reduced spectrum: R{E) is larger at low energies. For negative e. 
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R{E) increases with E and the slope is positive. In the small Am^ part of the LOW region 
the distortion of the spectrum is induced by the effect of the non-adiabatic edge (Fig. 1, 
Am^/4E = (0.3 — 3) x 10~^^ eV. Here the situation is opposite to that for the HIGH region. 
The slope is positive for positive e and negative for negative e as illustrated in Fig. 12(b). 

As follows from Figs. 11 and 12, the measurement of the electron energy spectrum 
provides information mainly on the value of Am^ but it will be difficult to determine e 
by measuring the slope (first moment of the spectrum) in the interval —0.3 < e < 0.3 
at Super-Kamiokande. At present Super-Kamiokande has presented the measured value 
(T) = 8.14 ± 0.02 MeV. The precision of this measurement is to be compared with the 
maximum theoretically expected variation A(T) = |(r)(e = 0.3) - {T){e = -0.3) | < 0.025 
MeV which occurs for the two values of Am^ shown in Fig. 12. Thus with the existing 
experimental precision, in the MSW region, the full range of e is allowed at ~ la. 

V. TESTS OF MAXIMAL MIXING IN FUTURE EXPERIMENTS 

In this section we consider the prospects of testing (near-)maximal mixing of the electron 
neutrinos in future experiments. We study various possibilities to measure the deviation e 
and we estimate the accuracy at which maximal mixing can be established or excluded. 
The extent to which the results of this Section hold in a three generation framework with 
\Ue3\ 7^ is discussed in Sec. VI. 

A. General requirements 

There are two requirements for a precise determination of the mixing: 

1. Uncertainties in the original neutrino fluxes should not play a role. For this purpose 
we will consider SSM independent observables, or at least observables which do not depend 
on the uncertainties in the boron neutrino flux. 

2. At least two independent observables should be measured. As we have seen in Sec. HI, 
the survival probability Pee and consequently all observables depend on both e and Am^. 
Even in the case of maximal mixing when the solar survival probability is constant, Pd — 1/2, 
a dependence of Pee on Am^ appears due to the Earth regeneration effect. 

Thus, to determine the mixing, one should find two independent observables which (i) 
are free of flux uncertainties, (ii) can be measured with high accuracy, and (iii) depend on 
different combinations of e and Am^. In what follows we will identify such observables and 
study the accuracy at which mixing can be measured. 

B. GNO and Super-Kamiokande 

The main objective of the GNO experiment [31] is to reach an accuracy ~ 3 SNU in 
the measurement of the Ge-production rate, Qoe- Also seasonal variations of Qce will be 
studied. Super-Kamiokande will continue to collect data for at least 10 years. With an 
energy threshold as low as 5 MeV the accuracy in measuring the Day-Night asymmetry will 
improve to ~ 0.010 - 0.015. 
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Notice that, in the MSW region, An-d is mostly sensitive to Am^, whereas Qoe h^s 
strong dependence on e. Therefore the pair of observables (^N-DiQce) is, in principle, very 
useful for measurements of the oscillation parameters in the matter conversion region. 

In Fig. 13 we show simultaneously contours of constant An-d at Super-Kamiokande site 
and Qoe in the Am^ — e plane. The iso-contours for Qoe have been obtained for central 
values of the solar fluxes according to BP98. The theoretical (la) uncertainty is about ±2 
SNU. 

The iso-contours of Qq^ and An-d are nearly perpendicular to each other, which re- 
flects the fact that these observables depend on different combinations of the oscillation 
parameters. However, the accuracy of the present experimental results is not enough to put 
statistically significant bounds on the mixing. The present experimental la intervals are 
Qgc ~ 69 — 80 SNU and An-d ~ 0.01 — 0.06. The resulting constraints on the mixing param- 
eters can be read from Fig. 13: in the HIGH region, —0.2 ^ e and Am^ ~ (1.5—8) x 10~^ eV^, 
and in the LOW region -0.10 ^ e ^ +0.15 and Am^ ~ (0.5 - 2) x 10"'^ eVl Inclusion of 
the theoretical uncertainties will expand these regions, especially in the HIGH case where 
the sensitivity of Qce to e is rather low. 

Notice that the present Gallium data somewhat disfavor maximal mixing in the HIGH 
region. We estimate that the la accuracy of the determination of e is 

Ae 0.15 - 0.20, (la). (73) 

Within 2a experimental errors the allowed regions cover most of HIGH parameter space of 
Fig. 13(a), and practically the entire LOW parameter space of Fig. 13(b). All values of 
e e (—0.3, -1-0.3) become allowed at the 2a level. 

With more data from GNO and higher statistics Super-Kamiokande measurements of 
the asymmetry one can reach better sensitivity. 

C. SNO 

The SNO experiment [36] will study various observables in three types of processes: 

1. Charged current neutrino-dcutcron scattering: the total rate above a certain threshold 
(we denote the reduced rate of events by [CC]), the energy distribution of events 
(electron energy spectrum) , and the time variation of events which can be characterized 
by the Day- Night asymmetry ^n-d; the zenith angle distribution, and the seasonal 
asymmetry. 

2. Neutrino-electron scattering: the total rate [ES], electron energy spectrum, and time 
variations. 

3. Neutral current neutrino-deuteron scattering: the total rate [NC], and time variations. 

Since the SNO observables depend on the boron neutrino flux only (we neglect the effect of 
the hep-neutrino flux), the ratios of rates are flux independent. Also relative time variations 
and energy spectrum distortion are flux independent. In what follows all the results will be 
presented for an energy threshold of 5 MeV. 
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In Fig. 14(a) we show the iso-contours of the double ratio [ES]/[CC] in the Am^ — e 
plane. In terms of averaged survival probability, P(Am^, e), the ratio can be written as 

[CC] P /(1-e) ^ ' 

The last equality apphes in the MSW region. The factor / = 0{1) and deviates from unity 
due to Earth matter effects which are small (see Eq. (34)) and r is the ratio between the 
— e and — e (ES) cross sections. The ratio [ES]/[CC] is larger than 1 and depends 
rather weakly on the oscillation parameters. From Eq. (74) we find the relation between 
the accuracy of measurements of [ES]/[CC], A([ES]/[CC]) and the corresponding accuracy 
of determination of e: 

Ae^^AfMV (75) 

That is, the accuracy is lowered by factor l/2r ^ 3. As follows from Fig. 14(b), in the 
MSW region the ratio [ES]/[CC] increases with e for fixed Am^. It varies within the limits 
[ES]/[CC]~ 1.15 ± 0.10 for e = 0.0 ± 0.3. This variation is comparable with the expected 
l(j error which is dominated by the uncertainty in the neutrino-deuteron cross-section (~ 
6%) [37,38] (statistical error has been calculated assuming 5000 CC and 500 ES events). 
Consequently no significant constraints on the oscillation parameters can be obtained, unless 
the uncertainty in the cross-section is reduced. According to Eq. (75), 10% precision in 
[ES]/[CC] corresponds to Ae k. 0.3. 

In Fig. 14(b) we show the iso-contours of the double ratio [NC]/[CC] in the Am^ — e 
plane. In terms of P(Am^, e), the ratio can be written as 

M.i = ^_ (76) 

icq p ■ ^ ' 

where, again, the last equality is valid in the MSW region. From Eq. (76) we get for the 
accuracy in the determination of e in the MSW region: 

Here the prefactor is smaller than 1. Moreover, the double ratio [NC]/[CC] will be deter- 
mined with much better accuracy than [ES]/[CC]. The theoretical uncertainties related to 
the neutrino-deuteron cross-sections essentially cancel out. The total la error, which in- 
cludes a statistical error for 5000 CC and 2000 NC events, is about 3.6%. According to 
Eq. (77), this corresponds to Ae ~ 0.04 for fixed Am^. As follows from Fig. 14(a), in the 
MSW region the ratio [NC]/[CC] increases with e for fixed Am^. It varies within the limits 
1.2 — 2.7 for e in the range (—0.3, +0.3). This variation is much larger than the expected la 
error, A([NC]/[CC]) ~ 0.08 (for [NC]/[CC] = 2). 

In the allowed regions of the oscillation parameters, the ratio [NC]/[CC] depends strongly 
on e. A precise determination of Am^ in these regions can be achieved from measurements 
of time variations, in particular, the Day- Night asymmetry ^n-d- ^ig. 15 we show iso- 
contours of A^Sj) and [NC]/[CC] in the Am^ — e plane. Notice that the asymmetry of the 
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CC-events is larger that the asymmetry of the ES-events for the same values of the oscillation 
parameters. The contours have weak dependence on e. The combined analysis of [NC]/[CC] 
(sensitive mostly to e) and ^n-d (sensitive to Am^) can give a precise determination of 
the oscillation parameters. According to this figure, measurements yielding [NC]/[CC]~ 
2 X (1 ± 0.04) and ^n-d ~ 0-1 x (1 ± 0.3) will determine e to an accuracy of order 

Ae K 0.05 - 0.07 (la) . (78) 

Notice that the same pairs of values (^n-d; [NC]/[CC]) appear in the HIGH and LOW 
regions. The LOW/HIGH ambiguity can be resolved by the KamLAND [39] reactor exper- 
iment which will give a positive oscillation signal in the case of the HIGH solution. It can 
also be resolved by the Borexino experiment which will show strong earth regeneration effect 
in the LOW region as discussed next. 



D. Borexino 

The Borexino collaboration [40] will measure the total rate of ES events and search for 
time variations of the signal. 

In Fig. 16 we show iso-contours of the reduced rate Rse (suppression factor) and iso- 
contours of the Day-Night asymmetry in the Am^ — e plane. As discussed before for Super- 
Kamiokande and SNO, the Day-Night asymmetry is sensitive mainly to Am^ while the 
deviation from maximal mixing can be restricted by the rate for which we can write 

Rbc = TBe + (1 - rBe)P , (79) 

where tbs ~ 0.24 is the ratio of the — e to z/g — e cross-sections for the beryllium neutrinos. 
From Eq. (79) and Fig. 16 we observe the following behaviours: 

(i) In the QVOl region, the survival probability P depends quadratically on e. We find 
Rbc ~ 0.62 with very weak e dependence. 

(ii) In the HIGH region, the transition between the quadratic and linear e dependence 
occurs. For Am^ ~ 10~^ eV^ the rate increases from 0.54 to 0.71 when e decreases from 
-1-0.3 to —0.3. The Day-Night asymmetry is very small here. 

(iii) In contrast, in the LOW region Borexino has higher sensitivity to the oscillation 
parameters. For maximal mixing, Rse is in the interval 0.65 - 0.70, and the asymmetry can 
be as large as 30%. The probabihty depends hnearly on e, so that we get from Eq. (79) 

Ae ~ ARbs. (80) 

1 - TBe 

Still, the sensitivity to e is low but Borexino will play a crucial role in fixing Am^ and 
in particular by the measurement of the Day-Night asymmetry will be able to resolve the 
HIGH/LOW ambiguity which may remain after the measurements at Super-Kamiokande 
and SNO. 
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E. LowNu Experiments 



A new generation of experiments aiming at a high precision real time measurement of 
the low energy solar neutrino spectrum is now under study [41]. Some of them such us 
HELLAZ, HERON and SUPER-MuNu [42] intend to detect the elastic scattering of the 
electron neutrinos with the electrons of a gas and measure the recoil electron energy and its 
direction of emission. The proposed experiment LENS plans to detect the electron neutrino 
via its absorption in a heavy nuclear target with the subsequent emission of an electron and 
a delayed gamma emission [43]. The expected rates at those experiments for the proposed 
detector sizes are of the order of ~ 1-10 pp neutrinos a day. Consequently, with a running 
time of two years they can reach a sensitivity of a few percent in the total neutrino rate 
at low energy, provided that they can achieve sufficient background rejection. This would 
allow the determination of e with a similar precision of a few percent, in particular in the 
QVOs, where SNO and Borexino cannot give much information due to their higher energy 
threshold. 

Sensitivity to the oscillation parameters in the LMA region can also be achieved in 
experiments detecting low energy z/g ffuxcs from nuclear reactors (detection threshold is 
about Etfi = 1.8 MeV). The Borexino experiment, in addition to detecting solar neutrinos, 
aims to detect the diffuse fluxes from nuclear reactors in Europe, mainly in Prance and 
Germany, at an average distance ~ 800 km. The KamLAND experiment [39] aims at 
detecting the low energy diffuse Pg fluxes from reactors in Japan from an average distance of 
~ 200 km. Both experiments can provide important information in discriminating between 
the HIGH and LOW solutions. However, they are expected to be very weakly sensitive to 
e. Given the short distances traveled by the neutrinos, matter effects are neghgible for both 
experiments. Consequently, the survival probability for these experiments takes the vacuum 
form of Eqs. (26) and (28) which depend only quadratically on e. With the expected 
achievable limits on the survival probability of about ~ 10-20%, a measurement of |e| < 0.3 
seems unfeasible. 

VI. EFFECTS OF THE THIRD NEUTRINO 

In the general case of three neutrino mixing when jC/esI 7^ the definition of deviation 
from maximal mixing becomes ambiguous. Formally, maximal mixing of the z/g in the three 
neutrino context can be defined as the equality 612 = vr/4 where 612 is the rotation angle 
in the plane of the first and second mass eigenstates. Phenomenologically, deviations from 
maximal mixing can be defined as a deviation of the averaged probability from 1/2 or a 
deviation of the depth of oscillations from 1. 

Let us consider neutrino mass spectrum which explains also the atmospheric neutrino 
problem via (mainly) — i/j- oscillations. This implies 

|Am|| = \Aml,\ < \Aml,\ ~ \Aml^\ = \AmlJ. (81) 

In this case (as far as solar neutrinos are concerned) the oscillations driven by Aml^^ can 
be averaged. This situation justifies our use of 775 as defined in Eq. (21) and tje as defined 
in Eq. (31) in the calculation of the survival probabilities. We get the following expression 
for the survival probability [44]: 
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Pee^{l-\Ues\YP2+\Ue3\\ (82) 

where P2 = P(Am^, 6, V) is the two neutrino mixing survival probabihty determined by 

Am^^Aml, UnO ^ \Ue2/Uei\, V ^ Vo{l - \U,s\'') . (83) 

where Vq is the matter potential Vq = \/2GFpYg/mN. In principle we can define a parameter 
62 which describes the deviation from maximal mixing of 9 defined in Eq. (83) similarly to 
the two neutrino parameters e and 6ex of Eq. (3). Then the deviation from maximal mixing 
in the three neutrino context will depend on the specific physical situation. If the mass 
splitting AttIq and the mixing angle 9 induce vacuum oscillations, then P2 ~ (1 + ei)/^ and 
Eq. (82) gives 

P,,^l + |-|[/,3|2. (84) 

Thus one can define the deviation from maximal mixing in the three neutrino context in 
this case as £3 = e| — 2|[/e3p- Clearly, the three neutrino case reduces to the two neutrino 
case if 

\Ue3\'^\el. (85) 

Taking |?7e3p at the level of the CHOOZ bound [16], |f/e3p ~ 0.05, Eq. (85) gives jes] > 0.3. 

If Am| and 9 lead to the adiabatic conversion in matter, then P2 ~ (1 — e2)/2 and 
Eq. (82) gives 

Pee~^-|-|C/e3r. (86) 

Here it is useful to characterize the deviation from maximal mixing through = e2+2|C/e3p. 
Corrections due to the third generation can be neglected provided that 

l^esP « ^£2. (87) 

For iC/eaP ~ 0.05, Eq. (87) gives €2 > 0.1. 

We have verified that, in the case of adiabatic transitions, the results of our calculations 
of the expected rates in the two-neutrino mixing scenario can be translated to a good ap- 
proximation to the case of three-neutrino mixing with the simple replacement of 62 with 
e^. This applies, for instance, to the contours for the Ar-production rate in Fig. 6 and the 
predictions in Fig. 7 and Fig. 8 in the range 10~^ > Am^/eV^ > 10~^. It also applies to the 
Gc-production rate in Fig. 13(b) and the corresponding predictions in Fig. 7 and Fig. 9(b) 
and 9(c). For the predictions of the SNO rates it can be used for Figs. 15 and 14 in the 
range 10^'' ^ Am^/eV^ ^ 10"'', and for the predicted rates at Borexino in Fig. 16(b) for 
Am7cV2 > 5 X 10-^. 

In the case of averaged vacuum oscillations (QVOl) the results for the three-neutrino 
scenario can be read from the results presented here with the replacement of 62 with eg. 
For long wavelength oscillations (QVOs and VO), the results for the three-neutrino scenario 
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cannot be directly derived from our results. Thus our predictions in these regions only hold 
for very small values of the mixing angle {Uesl, ^^^^ below the present CHOOZ bound. 

Concerning the value of Ues, although certain improvement on the present CHOOZ bound 
may be expected from long basehne experiments, such as K2K and MINOS [45], their final 
sensitivity is still unclear as it depends on their capability of discriminating against the i/e 
beam contamination. Ultimate sensitivity can be achieved at experiments performed with 
neutrino beams from muon-storage rings at the so-called neutrino factories [46]. 

VII. MAXIMAL MIXING AND OTHER EXPERIMENTS 

In the previous sections we have concentrated on effects in the solar neutrinos. Mixing 
of the electron neutrino can be probed in a number of other experiments. 



A. Atmospheric neutrinos 

Maximal and near-maximal mixing of the electron neutrinos can be probed in the at- 
mospheric neutrino studies. The oscillations in the Earth matter with parameters from 
the LMA or HIGH regions can give an observable effect in the e-like events. The electron 
neutrino flux at the detector [47-50] can be written as [50] 



Pe^{r COS^ 023 



(88) 



where r = F'^/F^ is the ratio of the original electron and muon neutrino fluxes, Pg/^ = 
P(AmQ, €2) is the two neutrino transition probability, and ^23 is the — Vt mixing re- 
sponsible for the dominant mode of the atmospheric neutrino oscillations. The transition 
probability can be of order one at Am^ > 3 x 10~^ eV^. It decreases fast with Am^/E" 
due to matter suppression of the mixing. Thus the biggest effect is expected in the low 
energy (sub-GeV) events sample. Notice that the probability Pg/* enters in Eq. (88) with 
a "screening factor" (rcos^^23 — 1) which turns out to be small. Indeed, for the sub-GeV 
sample r ~ 2 and the screening factor is exactly zero for maximal mixing in the atmospheric 
neutrinos [50,51]. The factor equals approximately cos 2^23, so that for sin^ 2^23 = 0.95 we 
get about 0.22. 

For ^^23 < 7r/4 the oscillations lead to the excess of the e-like events. Indeed some excess 
is hinted by the SK data. The excess can be deflned as Ne/N^ — 1, where A^'e and are 
the numbers of events with and without oscillations. This excess can be written in a matter 
dominant regime of oscillations {tje ^ 1) as 

^ - 1 R:; cos 2^23^^1(1 - 62 Ve), (89) 

e 

where 62 was deflned below Eq. (83). The excess depends on 62 linearly and it increases with 
r]E- However it is even more sensitive to deviations of ^^23 from the maximal mixing value. 
For Aml-^ = 10"^ eV^ and sin^ 2^23 = 0.95 the excess reaches about 3% and the dependence 
on e is weak. For Amli = 2 • 10""^ eV^ and the same 6'23, we flnd the excess about 4.5 %. 
Signiflcant dependence on e appears for Amfa = 3 x 10~^ eV^. However, it is unhkely given 
the size of the effect that atmospheric neutrino data will give any signiflcant information on 
the value of €2- 
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B. Supernova neutrinos 



Maximal or near-maximal mixing of the electron neutrinos will significantly modify 
properties of neutrino bursts from supernovae. The effects depend crucially on features of 
the whole neutrino mass spectrum and in particular on the value of C/gs and whether the 
mass hierarchy is normal (Am|i > 0) or inverted (Am|]^ < 0). Let us summarize here the 
main results (for more details sec [52] and references therein). 

All oscillation and conversion effects in supernova neutrinos are determined by the total 
survival probability of the electron neutrinos which in this subsection we will write as p, and 
total survival probability of the electron antineutrinos, p. (This property is related to the 
fact that the original spectra of the muon and the tau ncTitrinos arc identical and that the 
muon and tau neutrinos cannot be distinguished at the detection point). 

The probabilities should include the effects of propagation inside the star, on the way to 
the Earth and inside the Earth. Using the probabilities p and p, one can write the fluxes 
of the electron neutrinos, Fg, and electron antineutrinos, Fg, at the detector in terms of 
the original electron ( ant i) neutrino fluxes, F^ and Fg, and the non-electron neutrino flux 

pO = pO ^ PO ^ PO ^ po. 

X — n T IX T ' 

F,=p-F^ + {l-p)-F^ , F, = p-F^+{l-p)-F^ . (90) 

In general, p, and p depend on the neutrino energy. 

Let us summarize the results for specific neutrino mass and fiavor spectra. 

1) If |^7e3p > 3 X 10~^ the conversion in the resonance related to the largest (atmospheric) 
splitting (Am|i) will be completely adiabatic and the final effect depends on the type of 
mass hierarchy. In the case of normal hierarchy {u^ is the heaviest state) the resonance 
conversion occurs in the neutrino channel and for the survival probability we get 

P^|f/e3p<L (91) 

This probability is practically independent of the properties (mass, ffavour) of the first and 
second mass eigenstates. In particular, there is no sensitivity to e and no Earth matter effect 
is expected for neutrinos. 

In contrast, the antineutrino channels will not be affected by the high resonance and p 
will be determined by physics of the two light levels. 

For parameters in the HIGH and LOW regions, the neutrino propagation in the star is 
adiabatic, so that the survival probability in the star equals 

p^cos^9i2^\{l + e). (92) 

This probability can be further modified due to oscillations in the matter of the Earth. Thus, 
we expect the following consequences: (i) disappearance of the neutronization peak; (ii) 
hard v^. spectrum (coinciding with the original v^) spectrum at the cooling stage: 

Fe ~ Fl ■ (93) 

(iii) composite v^. spectrum: 
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F,= l{F^ + F^)-'-{F'^-F^y, (94) 

(iv) strong Earth matter effect (wfiidi leads to different signals at various detectors). For 
the HIGH mass range, the Earth effect is maximal in the high energy part of the spectrum, 
£^ > 20 MeV, whereas for the LOW solution the largest effect is in the low energy part. 

According to Eq. (94), the e-dependent term is proportional to the difTerence of the 
original fluxes. Thus due to the uncertainties in the predicticted fluxes it will be difficult 
to measure e. In order to reduce the theoretical uncertainty, one could in principle compare 
numbers of z/g and z/g events at large E which are determined by, respectively, F° and 
F°(l — e)/2 and are proportional to the same flux. 

For parameters of the two light states in the VO region, the neutrino propagation in the 
star is non-adiabatic, so that the survival probability can be writen as 

p^{l-Pc) cos^ 9i2 + Pc sin^ = ^ + e Q - P^) , (95) 

where the jump probability Pc depends on the details of the density profile in the outer 
regions of the star and cannot be rehably predicted. Practically, the probabihty should he 
between the adiabatic value (92) and the pure vacuum oscillation expression p — {1 + e^)/2. 

In the case of inverted mass hierarchy the sensitivity to e appears in the neutrino chan- 
nel and neutrinos and antineutrinos interchange their roles. Now the resonance is in the 
antineutrino channel so that 

P~|C/e3|'<l, (96) 

and the oscillations in the neutrino channels will be determined by physics of the two light 
levels. For parameters from the HIGH and LOW regions the propagation in the star is 
adiabatic and 

pRisin2^i2 = ^(l + e). (97) 

Moreover, Earth matter effects are expected for neutrinos. For the VO region we find, 
similarly to Eq. (95), 

P^l-'(l-Pc). (98) 



2 V2 



In this inverted scheme we predict: (i) partial disappearance of the i/g neutronization 
peak; (ii) hard spectrum of the electron antineutrinos: 



F, ^ ; (99) 



(iii) composite spectrum: 



Fe-l{F^ + FS)+'-{F'^-F^); (100) 



(iv) Earth matter effects are expected in the neutrino channel only. 
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2) If |f/e3p < 3 X 10^®, the effect of the third neutrino can be neglected: in the resonance 
channel the transition driven by Ue3 is strongly non-adiabatic, and in the non-resonance 
channel the mixing is always very small. In this case the problem is reduced to transitions 
in two level system with parameters determined by Am^i = Am^ and e. As a result both 
neutrino and antineutrino channels turn out to be sensitive to e. The effects include those 
considered above both for normal and inverted mass hierarchy. 

For HIGH and LOW regions of parameters the propagation proceeds adiabaticaly, and 
for the survival probabilities we get the expressions given in (92) and (97). Correspondingly, 
neutrino and antineutrino spectra will be given by Eq. (94) and (100). Thus we predict 
that both neutrino and antineutrino spectra will be composite, consisting of nearly equal 
admixture of the soft and hard components. In the high energy part where effects of the 
soft components can be neglected we get from Eqs. (94) and (100): 

§~l + 2e. (101) 

That is, larger signal (as compared with z7g) is expected for e > and smaller for e < 0. 
Also Earth matter effect is expected in both neutrino and antineutrino channels. 

3) If \Ue3\^ is in the intermediate region, 3 x 10~^ — 3 x 10~^, the adiabaticity in high 
mass resonance is partially violated and we expect some intermediate situation between those 
described in 1) and 2). In particular, both and z/g spectra will be composite, however 
admixtures of the soft and hard components will be unequal, etc.. 

To conclude, one expects strong influence of maximal and near-maximal mixing on the 
properties of the neutrino bursts. However, the uncertainties in the predicted neutrino 
spectra will make it difficult to obtain high sensitivity to e. 

Notice also that the analysis of the SN1987A data gives the 99% CL bound on p > 0.65. 
This would correspond to p > 0.3 [53]. Some recent calculations show that the difference 
between Pg and original spectra can be rather small, which would somewhat relax the 
above bound. 



C. Neutrinoless double beta decay 



The effective Major ana mass of the electron neutrino rriee relevant for the neutrinoless 
double beta decay is sensitive to the distribution of the electron neutrino flavor in the mass 
eigenstates (see [54] for recent discussion). The contribution to the effective mass niee from 
the two mass eigenstates responsible for the solar neutrino conversion can be written in 
terms of the oscillation parameters as 



mi(l + e) + Jml + AmKl — e)e 



(102) 



where mi is the mass of the first eigenstate and 0i2 is the relative phase of the first and 
tlie second mass eigenvalues. In the case of strong mass hierarchy, ml <S Amg we get from 
Eq. (102) 



(103) 
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According to this equation in the HIGH region the effective mass can be as big as (1—2) x 10~^ 
eV which can be probed at the next generation of the double beta decay experiments [55]. 
Notice that the contribution from the third mass eigenstate is strongly restricted by present 
experimental bound rUee < 2 x 10~^ eV. 

Although the dependence of rriee on e is rather strong, it will be difficult to measure e 
due to the uncertainties in the nuclear matrix elements. Eq. (103) can be considered as a 
test equation: if the measured values of mge, e and AttIq indeed satisfy this equation (within 
experimental and theoretical uncertainties) it will testify for the validity of whole scheme. 

In the case of strong mass degeneracy, m\ ^ Am^ we get 



1 

mee ~ -'Til 



1 + e) + (1 - 6)6'"^^^ , (104) 



where both non-oscillation parameters mi and </)i2 are unknown. For </)i2 = vr and (0) the 
mass equals rriee = triie (mi), so that e determines the lower bound on rriee- If nT-ee and e 
are measured, the above equalities will determine the upper bound and the lower bounds on 
the absolute scale of the neutrinos mass: rriee < rrii < niee/ e- 

In the case of inverted mass hierarchy the two states responsible for the solar neutrino 

conversion are degenerate: mi ~ m2 ~ \J Amatm ^-nd the effective majorana mass can be 
written as 



^ " 1 + e) + (1 - e)e'^^^ . (105) 



2 

In this case the measurement of e will allow us to determine the phase 0i2. According to 
Eq. (105), ^Aml,^e < niee < \J Am^tm which can be used as a test inequality for a given 
scheme. 

Thus measurements of e in the oscillation experiments will allow to determine or restrict 
the effective mass mge in the context of certain schemes of neutrino masses and mixing. 



VIII. CONCLUSIONS 

In this work we have explored the phenomcnological consequences of (near-) maximal 
mixing of electron neutrinos with other standard neutrinos. The possibility of such maximal 
or near-maximal lepton mixing constitutes an intriguing challenge for fundamental theories 
of flavour. Our aim was twofold. First we have formulated the present status of maximal 
mixing of i/^ in the light of existing experimental data from solar neutrino experiments. 
Second we have explored the best ways to measure deviations from such maximal mixing at 
future experiments. 

We show in Sec. II that both probabilities and observablcs depend on e quadratically in 
the regions of Am^ where the effects are due to vacuum oscillations, and they depend on e 
linearly when matter effects dominate. Consequently, for |e| <^ 1 the highest sensitivity to 
deviation from maximal mixing can be achieved in the Am^ ranges of the MSW effect. 

The results of a global fit to the existing solar neutrino data are presented in Sec. HI and 
summarized in Figs. 2 —5. From this analysis we find that values of the mixing parameter 
|e| = |1 - 2sin^6'| < 0.3 are allowed at 99% or lower CL for Am^ > 1.5 x 10"^ eV^ (which 
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contains the HIGH and QVOl regions) and for 4 x 10"^° eV^ ^ Am^ ^ 2 x 10"^ eV^ (which 
contains the defined LOW, QVOs and upper VACl regions). 

The role of the individual existing experiments on the determination of these regions is 
discussed in Sec. IV. We conclude that the present sensitivity to the mixing angle arises 
from the measurements of total event rates. The present data from Homestake experiment 
in Ar-production rate gives the strongest constraint on maximal or near-maximal mixing as 
it favours a significant deviation from e = 0. This conclusion is independent of the existing 
theoretical uncertainty on the boron fiux as discussed in Sec. IV B. On the other hand the 
measurement of the zenith angle dependence and the recoil electron energy spectrum are 
important in the determination of the allowed mass ranges but are very weakly sensitive to 
deviation from maximal mixing. With the present existing sensitivity all values of |e| < 0.3 
are allowed within 4<j. 

In Sec. V we have discussed the ways to improve our knowledge on deviations from 
maximal mixing at future experiments. We concentrate on observables which are SSM (or 
at least boron flux) independent. In order to determine both the mass and the mixing we 
study pairs of observables. First we have looked at the maximal sensitivity which may be 
achievable on the presently running experiments GNO and Super-Kamiokande. In principle 
the measurement of the Ge-production rate at GNO and the Day-Night asymmetry at 
Super-Kamiokande and SNO can give crossed information on the oscillation parameters in 
the matter conversion region. In practice, however, the expected sensitivity is not enough to 
substantially improve the present knowledge on Am^ and e. The role of SNO and Borexino 
experiments is discussed in Sees. VC and VD. We show that with the expected theoretical 
and statistical uncertainty the most sensitive observable to the mixing angle is the rate 
[NC]/[CC] measurable at SNO. For instance, a measurement yielding [NC]/[CC] ~ 2 x (1 ± 
0.04) and ^n-d ~ 0.1 x (1 ± 0.3) wiU determine e to an accuracy of order Ae ~ 0.07. There 
exist however an ambiguity on the allowed mass range between HIGH and LOW regions. We 
show that the LOW/HIGH ambiguity can be resolved by the measurement of the Day-Night 
asymmetry at Borexino experiment which is sensitive to strong Earth regeneration effect in 
the LOW region or by the detection of oscillations in long baseline reactor experiments such 
as KamLand. However no substantial improvement on the knowledge of e is expected neither 
from Borexino nor from the new generation of low energy experiments either with solar or 
reactor neutrinos. 
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Anr (eV') 

FIG. 5. The dependence of Ax^ on Am^ for e = —0.3 (dashed line), (solid line) and 0.3 
(dash-dotted line). Dotted horizontal lines mark 99% CL. Below this line the corresponding oscil- 
lation parameters are accepted at a confidence level lower than 99%. The three panels correspond 
to the following fits: (a) only total rates; (b) all the data except for the Homestake result; (c) all 
data with the boron neutrino flux treated as a free parameter. 
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Am' (eV') 

FIG. 7. (a) The Germanium production rate, (b) the Argon production rate and (c) the rate of 
events at Super-Kamiokande, as functions of Am^ for 3 values of e: +0.3 (dash-dotted hues), 
(soHd hues), —0.3 (dashed hues). The rates are normahzed to the BP98 no oscillation expectation. 
Also shown are the experimental constraints and the SMA prediction. See the text for more details. 
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FIG. 8. The dependence of the Ar-production rate on e for various values of Am^ within the 
(a) HIGH, (b) LOW, (c) QVOs and (d) VO regions. Full, dashed and dotted lines correspond to: 
(a) Am^ = 4, 0.8, and 0.2 x 10"^ eV^, (b) Am^ = 1, 0.7, and 0.3 x lO"'^ eV^, (c) Am^ = 10, 5, 
and 1 X 10~^ eV^, (d) Am^ = 5, 0.3, and 1 x 10^^° eV^, respectively. For each mass the three 
curves correspond to the predicted central value plus and minus la theoretical uncertainties. The 
boron neutrino flux is normalized to the measured rate at Super-Kamiokande. The horizontal 
band corresponds to the Homestake experimental result and its la error. 
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FIG. 9. The dependence of the Ge~production rate on e for various values of Am^ within the 
(a) HIGH, (b) LOW, (c) QVOs and (d) VO regions. Full, dashed and dotted lines correspond to: 
(a) Am^ = 4, 0.8, and 0.2 x lO'^ eV^, (b) Am^ = 1, 0.7, and 0.3 x lO"'^ eV^, (c) Am^ = 10, 5, 
and 1 X 10~^ eV^, (d) Am^ = 5, .3, and 1 x 10^^'^ eV^, respectively. For each mass the three 
curves correspond to the predicted central value plus and minus \a theoretical uncertainties. The 
boron neutrino flux is normalized to the measured rate at Super-Kamiokande. The horizontal 
band corresponds to the averaged Gallium result and its \a error. 
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